Dynamic Scaling in Diluted Systems Phase Transitions: Deactivation trough Thermal 
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Activated scaling is confimied to hold in transverse field induced phase transitions of randomly 
diluted Ising systems. Quantum Monte Carlo calculations have been made not just at the percolation 
threshold but well bellow and above it including the Griffiths-McCoy phase. A novel deactivation 
phenomena in the Griffiths- McCoy phase is observed using a thermal (in contrast to random) dilution 
of the system 



The presence of quenched disorder in quantum phase 
transitions at zero temperature is a topic of current in- 
terest. Two of the most important special properties of 
the disordered quantum phase transitions are the appear- 
ance of activated dynamic scaling and the existence 
of GrifRths-McCoy singularities [nU3| even away from 
the critical point. 

Activated dynamic scaling was first analytically proved 
to hold in the disordered one dimensional Ising model in 
a transverse field [pj and there are many other results for 
this model in the literature [Hj|JoMp|. The disordered 
Ising model in a transverse field considered as a quan- 
tum spin glass provides a reasonable description of the 
system LiH OxYi^xFi and also, considering the exis- 
tence of long-range correlated disorder in it, may be an 
appropriate model to describe non-Fermi liquid behavior 
in certain /-electron systems [^,|l|,|2|. Activated dy- 
namic scaling implies the existence of an infinite dynamic 
critical exponent {z — oo). It means that instead of the 
typical power law relationship between the characteristic 
time scale and the characteristic length scale, ^r '^ C^, a 
new exponential relation appears ^r ^ exp{const x ^'^) 
with t/j — 1/2 for the one-dimensional Ising model 0. 
There are nearly no analytical results for higher dimen- 
sions. The activated dynamic scaling seems to disappear 
for two-dimensional p3[ and three dimensional Ising glass 



systems 14 1 , but it has been proved to hold in the disor- 
dered two-dimensional Ising model in a transverse field by 
means of Quantum Monte Carlo simulations |lq,&ql and 
renormalization group analysis p7[ |. The only analytical 
prediction of activated scaling in dimensions higher than 
one has been made for the dilution probability transition 
at the percolation threshold {pc) of a diluted Ising sys- 
tem in a transverse field ||l^, where percolation critical 
exponents have been found. 

The phase boundary of this model at T = was stud- 
ied a long time ago ||l^,g^,^ and is expected to have a 
multi-critical point and a straight vertical phase bound- 



ary separating the ferromagnetic phase from the quan- 
tum Griffiths phase at the percolation threshold. The 
existence of this boundary and the activated scaling pre- 
dicted by Senthil et al. has been recently checked by 
means of Quantum Monte Carlo simulations, comparing 
the dilution probability transitions at constant transverse 
fields for p values bellow an above the one corresponding 
to the multicritical point |23]. 

However, there is no study of the existence of acti- 
vated scaling for transitions tuned by the transverse field 
in diluted Ising systems. These systems present a great 
advantage: fixing the probability of occupied places (p) 
and varying de transverse field (F) at T = 0, it is pos- 
sible to follow the evolution of the activated scaling at 
zero temperature in the ferromagnetic phase (p > pc) at 
the phase boundary (p = pc) and in the Griffiths phase 
(p < Pc)- In this paper we will use the diluted Ising model 
in a transverse field to directly determine the existence 
and evolution of activated scaling in the three different 
zones of the phase diagram. We will show how the acti- 
vated scaling exists not just at the percolation threshold 
but well above it and how it reaches a maximum near the 
phase boundary {p = Pc) keeping mostly constant all the 
way through the Griffiths zone. We start from the two- 
dimensional (d = 2) diluted Ising model in a transverse 
field with Hamiltonian given by 
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where et £ {0, 1} are quenched random variables produc- 
ing the dilution. Si — 1 with probability p and e^ = 
with probability 1 — p. In order to perform the Monte 
Carlo calculations we use the Suzuki- Trotter decomposi- 
tion [ p3[ and we map the phase transition at T = of the 
two-dimensional quantum system in a three dimensional 
classical system with action 
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FIG. 2. {LT)m vs. L for different values of the spin con- 
centration. Thin straight lines are just a guide for the eye. 
Thick line represents the behavior expected for the pure sys- 
tem (z = 1). 



FIG. 3. ln{Lr)m vs. L for different values of the spin con- 
centration. Thin hnes are the hnear fittings and the thick Une 
is the behavior expected for the pure system (z = 1). Inset 
shows the behavior of the exponent V(p)- Dotted line is a 
sigmoidal fitting. 

So, basically, the results presented indicate that the 
scaling is activated not just near the percolation thresh- 
old but well above it, and that it keeps nearly constant 
when the transition considered is inside of the GrifRths- 
McCoy quantum zone, at least for the values of p con- 
sidered, which are near pc (it may be also a maximum). 
Note that we are using, up to now, a purely random 
procedure to produce the diluted samples. 

The existence of the Grifhths-McCoy zone is due to 
rare regions which are locally in the wrong phase. Basi- 
cally there are some strongly coupled regions or clusters 
in the ferromagnetic phase even for spin concentrations 
smaller than the percolation threshold of the system. 
This zones are found studying the tail of local suscep- 
tibilities |0,||,|9|,|o|j3l]j3|,^. The quantum transitions 
allow the existence of these clusters in the wrong phase 
but they are very difhcult to detect in classical systems. 
As we have shown the phase transitions (tuned by the 
transverse field) due to these rare regions also present 
activated scaling. 

However there might be a possible way to deactivate 
the scaling in the Griffiths zone. Using the Suzuki- 
Trotter formalism the existence of the wrong phase comes 
from clusters which are infinite in the imaginary time 
direction, but they still are finite in the perpendicular 
slices. That is the reason why is possible to "feel" the 
effect of the dilution by the existence of an activated scal- 
ing. However, if the dilution is not introduced randomly, 
but incorporating an infinite correlation lengtii be- 



tween spins and vacancies some of the clusters formed 
in the slices may be considered to have infinite size. 
The wrong phase arising in these clusters is equivalent 
to the pure system ordered phase and is not affected by 
disorder. In this case the system in the Griffiths zone is 
expected to show deactivated scaling. Of course, this will 
happen just if the phase transition is due to the existence 
of strongly coupled regions (i.e. if the system is in the 
Griffiths phase {p < Pcj), but if the magnetic response is 
due to whole system (i.e. p > pc) the introduction of long 
range correlated disorder will affect changing the univer- 
sality class of the system (as in thermal classic transitions 
psf ) and enhancing the activation of the scaling [|2| . 

A way to produce this kind of dilution is using the 
tiiermal dilution instead of a random dilution Q] . With 
this kind of dilution the system belongs to the universal- 
ity class of long-range correlated disordered sys- 
tems with an exponent a = 2 — 77 |35|| , being 77 equal 
to 0.25 for the classical [d = 2) Ising model (the way 
to produce thermal dilution is described in |3J]). Basi- 
cally the procedure is as follows, the pure system is first 
thermalized to criticality and then one kind of spins are 
turned into vacancies. The samples produced by ther- 
mal dilution at criticality will have an spin concentration 
near p = 0.5, so we can compare them only with samples 
diluted randomly with probability p = 0.5. Both kinds of 
dilution will be at the Griffiths phase since in both cases 
p,(a)>0.5|p. 

The analysis performed has been exactly the same as 
before but now going up to L = 40. By the Binder 
Cumulant we have been able to determine the critical 
transverse magnetic field (Tcithermal) ~ 3.13) and the 
relation between {Lr)m and L. Fig. 4 compares de be- 
havior of the thermal dilution with the behavior of the 
random dilution with p = 0.5. The thick line repre- 
sents the pure behavior (excepting proportional factors). 
Clearly the scaling has been deactivated, as expected, 
and z ~ 1 as corresponds to the behavior of the pure 
system. This phenomenon could never happen with a 
system away from the quantum Griffiths zone. 

In conclusion. Quantum Monte Carlo calculations in 
diluted Ising models in a transverse field at T = show 
that dynamic scaling holds above the percolation thresh- 
old {p > Pc), at the percolation threshold {p = Pc) and 
inside the quantum Griffiths-McCoy phase (p < Pc)- The 
evolution of the activated scaling has been characterized, 
showing how it grows monotonously towards the value 
corresponding to the percolation threshold and how it 
appears to remain nearly constant in the Griffiths zone 
for values of p near pc- A new way to deactivate the 
scaling in the Griffiths zone by thermal dilution has been 
proposed. The expected deactivation works due to the 
fact that phase transitions come from strongly coupled 
regions, and not from the whole sample. The deactiva- 
tion has been clearly confirmed by means of Quantum 
Monte Carlo simulations. 



FIG. 4. (I/T)m vs. L for a randomly diluted system with 
spin concentration p = 0.5 and a thermally diluted system. 
The thin line is a guide for the eye, and the thick line rep- 
resents the behavior expected for the pure case {z = 1) ex- 
cepting proportional factors. Note how the scaling has been 
deactivated for the thermal case. 
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